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Abstract 

First-principles density-functional calculations are performed to investigate the thermal transport 



i-rt ' properties in graphene nanoribbons (GNRs). The dimensional crossover of thermal conductance 

C ■ 

Q ■ from one to two dimensions (2D) is clearly demonstrated with increasing ribbon width. The thermal 



conductance of GNRs in a few nanometer width already exhibits an approximate low-temperature 



^ I dependence of T , like that of 2D graphene sheet which is attributed to the quadratic nature of 

a^ ■ 

dispersion relation for the out-of-plane acoustic phonon modes. Using a zone-folding method, 

oo : 

CN . we heuristically derive the dimensional crossover of thermal conductance with the increase of 

(<— ^ I ribbon width. Combining our calculations with the experimental phonon mean-free path, some 

(N ■ 

typical values of thermal conductivity at room temperature are estimated for GNRs and for 2D 

graphene sheet, respectively. Our findings clarify the issue of low-temperature dependence of 

Xj, 

^ , thermal transport in GNRs and suggest a calibration range of thermal conductivity for experimental 

measurements in graphene-based materials. 



I. INTRODUCTION 

Since the first exfoliation of graphene,- many exotic properties have been discovered in 
the two-dimensional graphene sheet.- Besides graphene, the electronic states of graphene 
nanoribbons (GNR),-~- thin strips of graphene, can be tuned through controlling the ribbon 
width^"^ and the edge chirality.-"^ Graphene-based materials can be a technological alterna- 
tive to silicon semiconductors because of their unique properties.^"^ Apart from their elec- 
tronic properties, graphene turns out to be an excellent heat-conduction material.-"— How- 
ever, some discrepant values of thermal conductivity are experimentally reported. Thermal 
conductivities for the suspended single-layer graphene flakes are first measured^ using Ra- 
man spectroscopy with values in the range from 4.84 ±0.44 x 10^ to 5.30 ±0.48 x 10'^ W/mK 
at room temperature. Using a similar Raman scattering method or utilizing the conven- 
tional heat bath method, successive measurements on graphenes^"— have reported thermal 
conductivity varying from ~ 170 W/mK— to ~ 5000 W/mK-^i^. Monolayer of graphene 
in contact with silicon dioxide has a thermal conductivity of about 600 W/mK.- Room- 
temperature thermal conductivity^^ decreases from 2800 to 1300 W/mK as the number of 
atomic planes increases from two to four in few-layer graphenes. Some different theoretical 
results of thermal conductivity with various approaches have also appeared.—"— 

Furthermore, there is another confusing with regard to the low-temperature dependence 
of thermal transport in two-dimensional graphenes. For one- dimensional quantum atomic 
chains,— thermal conductance at low temperatures is proportional to temperature T with 
the quantized universal coefficient n'^k'^/Sh, where ks is the Boltzmann constant and h is 
the Planck constant. For two-dimensional materials, it is well known that thermal conduc- 
tance depends on temperature as Debye T^ law. A calculation by the valence-force field 
demonstrates that thermal conductivities at low temperature conform^^ to the conventional 
relation of T^. But a theoretical estimations^ and the continuum model^^ propose a T^'^ 
temperature dependence of thermal conductivities for graphene, followed by some recent 
experimental observations.—"— 

To understand thermal transport in graphene-based materials is not only important to 
the possible technological applications but also crucial for a fundamental understanding of 
thermal properties in low-dimensional systems. The reduced-dimension effects on thermal 
transport are becoming critical to both device reliability and intrinsic physics. Temperature 



dependence of thermal transport, especially in the low-temperature range, is essential for 
understanding the fundamental physics behind the measured thermal conductivities. An 
accurate knowledge of the temperature dependence of thermal transport is important be- 
cause it helps to identify phonon branches which have the dominant contributions to heat 
conduction. Low-frequency acoustic phonon branches which determines the behavior of 
low-temperature thermal transport are sensitive to the weak-bonding interactions between 
the long-range atoms. The valence-force field considering only the nearest neighbors or the 
continuum model may not be an appropriate model. Moreover, the experimental measure- 
ments of temperature dependence alone cannot provide evidence in favor of one or the other 
phonon contribution because the temperature dependence in graphite depend strongly on 
the material quality.^^ In contrast, the first-principles approach can yield accurate phonon 
dispersioii^^ relations without empirical parameters, including the long-range interactions 
among the atoms. 

In this paper, we report a systematic investigation on temperature dependence of ther- 
mal transport in graphenes varying from quasi one-dimensional nanoribbons to infinite two- 
dimensional sheets via first-principles density-functional theory calculations.— The dimen- 
sional crossover of thermal conductance from one to two dimensions is clearly observed in 
GNRs with increasing ribbon width. We find that the thermal conductance of GNRs of a 
few nanometers wide already exhibits an approximate low-temperature dependence of T^'^. 
Ballistic thermal conductance of wide ribbons is also found to converge towards the cor- 
responding values of 2-D graphene sheet at low temperature. Combining our calculations 
with the experimental phonon mean-free path(MFP),- typical values of thermal conduc- 
tivity are estimated at room temperature for GNRs and 2D graphene sheet, respectively. 
We conclusively clarify the issue of low-temperature dependence of thermal conductance 
for graphene-based materials and suggest a calibration range of thermal conductivity for 
experimental measurements. We expect that our first-principles results may provide some 
valuable insights into the experimental measurement of thermal transport in graphene-based 
materials. 



II. COMPUTATIONAL METHODS AND FORMULAS 

To calculate the thermal conductance, we should obtain force constants in graphene 
nanoribbons and in 2D graphene sheet. The nanoribbons are classified^ by the number of 
dimers in unit cell using a common convention, where armchair-edged nanoribbons or zigzag- 
edged nanoribbons with n dimers are referred to as n-AGNR or n-ZGNR. To eliminate 
the dangling bonds at the edge of ribbon, we have passivated the nanoribbon edges with 
hydrogen atoms. The interatomic force constants are obtained from a series of supercell 
calculations using the software packages SIESTA— based on the first-principles density- 
functional approach. The ribbons are periodic with a supercell of 5 unit cells along the ribbon 
axis for minimizing the periodic continuation. In other two directions, the periodic image 
of ribbon is separated from its nearest neighbors at least by 18 A to prevent interactions 
between them. We use a 7 x 7 x 1 supercell for calculations of the graphene sheet, with 
a space 20 A between periodic images of sheet. For DFT calculations, Troullier-Martins 
pseudopotentials are used for both C and H atoms. The valence electrons are described by a 
double — C basis set plus polarization orbits (DZP). The local-density approximation (LDA) 
with the exchange-correlation functional due to Ceperley- Alder is employed. To suppress 
the force fiuctuations, we use a rather fine Mesh cutoff of 400 Ry during calculations. 

We first perform geometrical optimization with a force tolerance 0.001 eV/A with the 
aforementioned parameters using a conjugated-gradient minimization method. Force con- 
stants then are derived from the optimized structures with the method of finite difference, 
where each atom is displaced from its equilibrium position by a distance of ±0.0212A. After 
forces are evaluated, a central finite difference with respect to the displacement is utilized 
to compute the fore constants. We have calculated force constants for the n-AGNRs, n- 
ZGNRs with n = 2, 3, ■ ■ ■ ,14 dimers in width and for the two-dimensional graphene sheet, 
respectively. 

Thermal conductivity k along one direction, e.g. the z— axis direction, can be expressed 
by the well-known Boltzmann-Peierls formula^ 

:—fvjj,-——V^L, (1) 

where p denotes the polarization of phonon, q wave vector, uj^ the angular frequency, Jb 
the Bose-Einstein distribution at temperature T. The phonon group velocity f^ along the 



transport direction z is given by v^^ = dup/dqz and Ip represents the length of phonon 
mean-free path (MFP). The number n = 1,2 and 3 denotes the system dimension. The 
integration over the wave vector is carried out within the first Brillouin zone. Graphene 
has a long phonon MFP, which is estimated to be about 775 nm near room temperature.- 
Therefore, at low temperatures, it is reasonably expected that phonon transport in GNRs 
or in 2-D graphene sheets will be dominantly ballistic or quasi-ballistic rather than diffusive. 
In regime of ballistic or quasi-ballistic transport, thermal conductivity is not well defined 
because thermal conductivity k, depends on the length of system. Assuming that the MFP 
Ip for ballistic or quasi-ballistic transport is independent of the polarization since it can be 
compared with the length of system in experiments, we thus rewrite Eq. ([T]) using the mean 
value theorem as 



Kz = o'l, (2a) 

T^^^P^^l- (2b) 

Here / is the average length of phonon MFP, p denotes the polarization and a is defined 
as thermal conductance per unit area. In particular, thermal conductance a in Eq. ([2]) for 
one- dimensional ballistic transport is reduced to the following expression 

1 v^ /■ rfq dfB .„v 

where S = h ■ w is the cross-section, h the thickness and w the width. In contrast with 
thermal conductivity, thermal conductance is uniquely defined in the regime of ballistic 
transport. Here we concentrate on the salient features of temperature dependence of thermal 
transport at low-temperature in both GNRs and graphenes such that nonlinear scattering 
characterized by MFP is first neglected. Phonon transmission is assumed to be ballistic 
such that thermal transport at low temperature is distinguished by thermal conductance 
a as defined in Eq. ([2]). This is a valid approximation for thermal transport in graphenes 
at low temperatures, where the length of MFP is approaching the system size. Nonlinear 
scattering is later included through the average length of phonon MFP in the end of paper. 



III. RESULTS AND DISCUSSION 



A. Phonon dispersion relations 



Since thermal conductance is intrinsically related to phonon dispersion relations, we first 
present the calculated phonon dispersions from first principles. Fig. [1] shows the calculated 
phonon dispersion for the 8-ZGNR and 2-D graphene sheet. It can be seen from the figure 
that the stretching modes due to the C-H bonding in GNRs are located about 3110 cm~^. 
We can also find that the highest frequency of longitudinal-optical(LO) mode at F point 
is about 1650 cm~^ for both GNRs and the two-dimensional graphene sheet. For GNRs, 
there are four acoustic phonon branches near F point, including an out-of plane mode(ZA), 
an in-plane transverse mode(TA), an in-plane longitudinal(LA) and a torsion branch that 
will disappear with the increasing width. In comparison, similar acoustic branches except 
the torsion one are obtained for the 2-D graphene sheet. As shown in Fig. [H the ZA mode 
shows a quadratic energy dispersion around the F-point as a consequence of the point-group 
symmetry of graphene^ while the TA and LA modes display a linear dispersion. The phonon 
dispersion curves both for the GNRs and 2-D graphene sheet are consistent with the previous 
reported results.— 

The phonon dispersion between the GNRs and the 2-D graphene sheet can be mutually 
transformed through the unfolding of nanoribbons's Brillouin zone to that of graphene.— 
The phonon dispersions of the 2-D graphene sheet comprise three acoustic branches and 
three optical branches as indicated in Fig. Wijo) because there are two carbon atoms in the 
unit cell of graphene. In GNRs, a group of six modes can be found equivalent to the six 
phonon modes of graphene, with respect to the phonon eigenvectors near F — point. Hence, 
all phonon modes of GNRs can be interpreted as these six fundamental modes and their 
overtones.— 

B. Thermal conductance of graphene nanoribbons 

Using the complete phonon dispersion relations, we have systematically calculated ther- 
mal conductance a per unit area of Eq. ([2]) for GNRs with increasing width and for the 
2-D graphene sheet, as shown in Fig. O During the calculation, the layer thickness h for 
both GNRs and the graphene sheet is taken a typical value h = 0.335 nm. The ribbon 



width w for ZGNRs and AGNRs is computed by the formula i^zgnr = \/3(r2 — l)ao/2 and 
w^AGNR = {P' — l)ao/2, respectively. Here n is the number of dimers in unit cell and a^ denotes 
the graphene lattice constant given by Oq = 0.246 nm. For simplicity, a very small deviation 
of ribbon width after the structure relaxation is not considered in the formula. Fig. [2] shows 
thermal conductance per unit area as a function of temperature for the n-ZGNRs with the 
increasing number of dimers n. The open circles in Fig. [2] represent thermal conductance 
along T — K direction, corresponding to the axis direction of ZGNR. For a clear display of 
low-temperature dependence of thermal conductance, a log-log plot of thermal conductance 
for T < 200K is further illustrated in Fig. [31 

1. Dimensional crossover of thermal conductance 

The most significant feature of thermal conductance is the scaling behavior in the low- 
temperature dependence for ZGNRs. To quantitatively show such temperature dependence, 
the curves of thermal conductance a fitted vs temperature as T'' below 200 K are also 
plotted in Fig. |3l We can find that the exponent /3 increases from /3 = 1.09 to /3 = 1.41, 
approaching to the value of the two-dimensional graphene sheet /3 = 1.58, when the number 
of dimers for ZGNRs grows from n = 2 to n = 14. This dimensional crossover of thermal 
conductance from one to two dimensions is clearly demonstrated in Fig. [3t^a). It is well 
known that thermal conductance at low temperatures is proportional to temperature T for 
one- dimensional quantum atomic chains^ whereas the low-temperature dependence of T^ 
for the two-dimensional phonon gas is described by the conventional Debye law. For 2D 
graphene sheet, a calculation by the valence-force field demonstrates that thermal conduc- 
tivity at low temperature conform^^ to the conventional relation of T^. But a theoretical 
estimation^^ propose a T^'^ temperature dependence of thermal conductance for an infinite 
2D graphene sheet. An analysis of vibrational modes^ using the continuum mechanics 
shows that GNRs with width w > 500 nm will show a low-temperature dependence of T^'^. 
In comparison with the approaches of the empirical valence-force and continuum mechan- 
ics, the density-functional theory can yield an accurate phonon dispersion relations without 
empirical parameters. Our first-principles calculations shown in Fig. |3]^a) clarify that the 
two-dimensional graphene sheet will display a low-temperature dependence of T^'^ and that 
dimensional crossover of thermal conductance from one to two dimensions will be expected 



for the GNRs of only a few nanometers in width. We conclude that thermal conductance 
for most samples of graphene flakes in experiments with a width beyond a few nanometers 
will be characterized by a low-temperature dependence of T^'^. Our calculated results are 
consistent with the experimental results,—""— where the width of graphene samples is on the 
order of micrometers. 

Next we elucidate the reasons for the dimensional crossover of thermal conductance in 
GNRs through a heuristic derivation. The phonon wave vector q = {qt,qi) for GNRs consists 
of the transverse direction component qt perpendicular to the ribbon axis and the longitu- 
dinal component qi along the ribbon axis. The longitudinal wave vector qi represented by 
the X-axis in Fig. [I](a) is continuous because the ribbon length goes to infinite. But the 
transverse wave vector qt is discrete, due to the finite ribbon width which only allows the 
standing waves with the boundary condition qtw = k ■ n, k = 0,1,- ■ ■ ,n — 1. Here w 
denotes the ribbon width and n is the number of dimers in width. Therefore, the transverse 
wave vector qt can only take the quantized values given hj qt = k- n/w. Each of the discrete 
value of qt corresponds to one branch of over-tone phonon mode.— The vibrational modes 
in GNRs shown in Fig. [TJ^a) can be classified^ into: six fundamental modes equivalent to 
the phonon modes in 2-D graphene sheet, over-tone modes of the fundamental modes and 
the C-H modes resulting from the passivation with hydrogen at edges. The contribution to 
thermal conductance from the C-H modes can be neglected at low temperature on account 
of their high frequencies. Hence, the summation p in Eq. ([3]) over different polarization 
branches can be separated into the terms of six basic modes plus their over-tone modes, ex- 
pressed by Sp — )■ Sp/Sfe=o,n-i- Here p' denotes the summation over six fundamental phonon 
modes. When the number of dimers n in unit cell is large enough, we can transform the 
discrete summation k into the integration over continuum qt, mathematically described by 
— Sfc=o,n-i ~^ f d^lf In other words, for the nanoribbon with enough width, thermal con- 
ductance for one-dimensional transport in Eq. (J3]) will transition into the two-dimensional 
expression a2D = j^p'7^a2 /'^Qt'^li^p' "ot^p' ^i^^^ ^^^ constraint of Vp, > 0, where h is the 
thickness. Here p' denotes the basic phonon mode, ranging from LA,TA,and ZA to LO,TO 
and ZO modes shown in Fig. W(Jo) for the 2-D graphene sheet. From our numerical calcula- 
tions, we find that the number of dimers n = 14 already leads to low-temperature thermal 
conductance fairly close to that of 2-D graphene, demonstrating a dimensional crossover 
from one to two dimensions. 
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2. Temperature dependence of thermal conductance 

Further we turn to the causes of the low-temperature dependence of T^'^ in 2-D graphene 
sheet through an analytic derivation. The key reason lies in the different contributions to 
thermal conductance made by each of the polarized phonon branch with different dispersion 
relations. Note that we only need to consider the acoustic branches during calculating 
thermal conductance at low temperature. It can be observed from Fig. [I](b) that both LA 
and TA branches exhibit the linear phonon dispersion relations at low-frequencies while the 
ZA branch display an approximately quadratic dispersion curve near F point. To illustrate 
the dispersion relations in full phase space, we have presented the three-dimensional view 
of the ZA branch of the first brillouin zone in Fig. El^b). A direct proof of contribution 
to thermal conductance proceed as follows. In the integration of Eq. ([2]), we can take the 
integration over frequency u in substitute of the wave vector q by introducing the delta 
function, 

a2D oc }_^J -^hwp^Vp, (4a) 

with T[uj] = ^ d'^c[6{uj - ujp)Vp. (4c) 

p ^ 

Here 6 is the dirac delta function. T[u}] is defined as the effect transmission, which means 
the number of phonon branches at a given frequency u for one-dimensional transport. For 
two-dimensional transport, it is straightforward to prove that T[uj] behaves as co^^"^ for 
quadratic phonon dispersion relation while T[u}] is proportional to u with regard to linear 
phonon dispersion relation. Assuming that T[uj] oc w", it can be easily derived that thermal 
conductance has low-temperature dependence as o"2D oc T^+". Therefore, the contribution 
to thermal conductance made by the ZA branch with quadratic phonon dispersion scales as 
T^-^ while thermal conductance resulting from the linear LA and TA branch increases with 
temperature as T^. 

The reason why whole thermal conductance at low temperature increases with temper- 
ature as T^'^ like that of the ZA branch roots in differences of the vibrational density of 
states(vDOS), which can be defined as vDOS[(X'] oc f d'^q6{u — Up). It can be easily verified 
that vDOS in two dimensions is proportional to u for the phonon branch with linear disper- 



sion while two-dimensional vDOS for the quadratic phonon branch remains to be constant. 
Accordingly, as the phonon frequency tends to zero, the number of vibrational modes ex- 
pressed by vDOS also decreases to zero for branches with the linear phonon dispersion. On 
the other hand, the number of vibrational modes for the quadratic phonon branch does not 
change with frequencies. At low temperature, only low-frequency phonons mainly contribute 
to thermal conductance. Therefore, the ZA phonon branch with the quadratic dispersion 
makes a dominant contribution to low-temperature thermal conductance due to its constant 
vibrational density of states. By contrast, the contributions resulting from the linear LA 
and TA phonon branches are trivial because of their diminishing vibrational modes at low 
frequencies. Namely, thermal conductance at low temperature for two-dimensional graphene 
sheet increases with temperature as T^'^ on account of the dominant contributions made by 
the ZA phonon branch. 

3. Convergence of low-temperature thermal conductance and anisotropy. 

We now discuss another interesting feature observed in Fig. |2l It can be seen from the fig- 
ure that thermal conductance of GNRs at low temperatures converges to that of 2D graphene 
sheet with the increase of ribbon width. A systematical investigation on thermal conduc- 
tance as a function of the ribbon width for the zigzag and armchair graphene nanoribbons 
at temperature T=200K is shown in Fig. HI We can find that thermal conductance per unit 
area both for ZGNRs and AGNRs decreases with the increase of ribbon width, converging 
to a stable value. At low temperatures, the value of convergence agrees well with that of 2D 
graphene sheet along the corresponding direction, as shown by the open circles in Fig. |2l 
Such agreement can be explained by the fact that with increasing width phonon dispersion 
relations of GNRs, especially at low frequencies, are approaching that of 2D graphene sheet 
as described in the foregoing paragraphs. By contrast, there are some value differences in 
thermal conductance at high temperatures between GNRs and 2D graphene. The devia- 
tion results from the scattering caused by the confinement from the edges. Because of the 
confinement of standing wave due to the ribbon edges, the phonon dispersion curves tend 
to be flat at high frequencies, yielding a small group velocity. Thus, the reduced group 
velocity for high-frequency phonon branches in GNRs leads to smaller thermal conductance 
at high temperature. In addition. Fig. H] shows that thermal conductance per unit area of 
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the ZGNRs is larger than that of AGNRs. This behavior originates from the anisotropy of 
phonon dispersion relationships in consistency with other simulation results.—"— 

C. Nonlinear scattering and thermal conductivity 

Finally, we include the nonlinear phonon scattering effects through the phonon MFP. 
The above calculations have not taken into account phonon-phonon scattering. At high 
temperatures, the nonlinear phonon scattering is inevitable for more realistic situations. 
However, to calculate nonlinear phonon scattering for GNR system with a few hundred 
atoms is a formidable task for the firs-principles approach and is also beyond the ability 
of any first-principles implementations. For simplicity, we phenomenologically introduce 
the phonon-phonon scattering into the present ballistic theory through the phonon MFP 
described by Eq. ([2]). It is estimated from experiments that the phonon MFP at room 
temperatur^ on average is about 775 nm for graphene samples. At T=300K, our calcu- 
lated thermal conductance per unit area for 2D graphene sheet is o"2d = 5.5 x 10^ W /m^K 
along r — i^ direction and the converged value of <Jzgnr = 3.4 x 10^ W/m^K is obtained 
for the ZGNRs as shown in Fig. |2j In comparison with the previous estimation of thermal 
conductance^'^ for 2D graphene at room temperature o"2_d = 3.1 x 10^ W/m^K given by 
the formula a = 0.6 x lO^T^/^W/m^K^/^, our results are qualitatively consistent with this 
estimation. Combining the calculated thermal conductance with the experimental phonon 
MFP, thermal conductivity at room temperature k,2D = 4263 W/mk for 2D graphene and 
i^GNR = 2635 W/mk for ZGNRs is achieved. Even though our calculated results cannot be 
directly compared with experimental results where some other factors exist, such as struc- 
ture defects, supporting substrate and surface tensions, etc., a calibration range of thermal 
conductivity in GNRs and two-dimensional graphene is suggested for experimental measure- 
ments from our calculations. We propose that the room-temperature thermal conductivity 
of GNRs may be estimated on the order of ~ 2600 W/mK for GNRs and of ~ 4200 W/mK 
for 2D graphene sheet if intrinsic phonon-phonon scattering is only considered. We think 
that the measurement of intrinsic thermal conductivity in GNRs or in 2D graphene sheet 
can be reduced to measuring the accurate phonon mean-free path.— 
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IV. SUMMARY AND CONCLUSIONS 

In summary, we have investigated low-temperature thermal transport in GNRs and 2D 
graphene sheet using first-principles density-functional theory approach. A dimensional 
crossover of thermal conductance from one dimensional GNRs to two dimensional graphene 
sheet is clearly demonstrated. We determine that thermal conductance of GNRs with a 
width of a few nanometers will exhibit a low-temperature dependence of T^'^. A heuristical 
derivation of thermal conductance by the zone-folding method is carried out to elucidate the 
reason for this dimensional crossover. The reason for the temperature dependence of T^'^ 
can be attributed to the quadratic nature of dispersion relation for the out-of-plane acoustic 
phonon branch. In addition, we find that low-temperature thermal conductance in GNRs 
converges towards the corresponding values of 2D graphene sheet with the increase of rib- 
bon width. Combining our calculations with the experimental phonon MFP, typical values 
of thermal conductivity for GNRs and for 2D graphene sheet are estimated at room tem- 
perature. Our findings conclusively clarify the issue of the low-temperature dependence of 
thermal conductance for GNRs and suggest a calibration value range of thermal conductivity 
in graphene-based material. We hope that our findings can offer some enlightening insights 
into the experimental measurement of thermal transport in graphene-based materials. 
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FIG. 1: (Color online). Plionon dispersion relations calculated using first-principles density- 
functional calculations, (a) The dispersion relation for an 8-ZGNR. The wave number q is given in 
units of the length of unit cell T . (c) The phonon dispersion curves plotted along high symmetry 
directions for a two-dimensional graphene sheet. 
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FIG. 2: (Color online). Temperature dependence of thermal conductance per unit area for ZGNRs 
with the increase of ribbon width n, which indicates the number of dimers per unit cell. The open 
circles in the figure represent thermal conductance along V — K direction, corresponding to the 
axis direction of ZGNR. 
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FIG. 3: (Color online), (a) The scaling behavior of thermal conductance with temperature for the 
2-ZGNR, 14-ZGNR and two-dimensional graphene sheet, respectively. The dashed, dotted, and 
solid lines in the figure are numerically fitted as T^ at low temperatures, (b) The 3D-view and 
its projection of phonon dispersion relations for the out-of-plane acoustic branch(ZA) in the first 
brillouin zone. 
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FIG. 4: (Color online). Thermal conductance per unit area as a function of the ribbon width for 
the zigzag and armchair graphene nanoribbons at temperature T = 2Q{)K. 
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